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In contrast to other atomic systems, in true muonium [pApT) the leading-order Z boson correc¬ 
tions to the hyperfine splitting are shown to be experimentally accessible in the near future. This 
contribution ( — 109 MHz) constitutes a necessary contribution to a full 0{mcA) calculation of the 
true muonium hyperfine splitting. This calculation would enable a number of possible to the muon 
problem to be constrained. Additionally, we compute the general expression for a pseudovector 
coupling to particle-antiparticle bound states at leading order, including the annihilation channel. 

PACS numbers: 36.10.Ee, 31.30.jr, 32.10.Fn, 24.80.-|-y 


I. INTRODUCTION 

The current state of flavor physics might be described 
as having a “muon problem”. Several muon observables 
have been shown to have varying levels of disagreement 
with Standard Model calculations. In Table U are listed 
a few of the most prominent signals that are in tension 
with expectations. 


Obs. 

Discrep. Ref. 

(5-2)^ 

2.9cr [Ij 

rp from p~p'^ 

7a [2j 

B+ K+p+p- 

2.6o- [3] 

h ^ pr 

2.5o- [4| 


TABLE I. Set of prominent muon signals that have disagree¬ 
ments with theoretical expectations. 

A strong candidate for shedding light on the muon 
problem is the bound state {p'^p~), dubbed “true 
muonium” Q . Simpler bound states like positronium 
(e+e”), hydrogen, and muonium (^+e“) have attracted 
significant attention as testing grounds for precision QED 
studies^, but are limited in their BSM discovery po¬ 
tential by either the small electron mass or large uncer¬ 
tainties from unknown nuclear structure effects. In con¬ 
trast, true muonium has a much larger reduced mass, and 
its QCD corrections are limited to the better-understood 
hadronic vacuum polarization effects, due to its leptonic 
nature. 

Unfortunately, true muonium has yet to be directly ob¬ 
served. The first reason is the technical difficulty of pro¬ 
ducing low-energy muon pairs, coupled with the bound 
state’s short lifetime (t Ri 1 ps), which presents an inter¬ 
esting challenge to experimenters. Numerous proposed 
methods of production channels have been discussed over 
the years in the literature [Mil , but until recently, none 
have been seriously developed. A second, more prosaic. 


* hlammiv@asu.edu 


reason true muonium has been neglected is that, until 
the advent of the (g — 2)^ anomaly, there was no ap¬ 
parent reason to expect that true muonium would offer 
any novel physics compared to the simpler, more easily 
produced bound states. 

Currently, the Heavy Photon Search (HPS) lr^ e xner- 
iment has plans to search for true muonium IlSlI. and 
DImeson Relativistic Atom Complex (DIRAC) [1^ has 
discussed the possibility of its observation in an upgraded 
run[l^. In both situations, the true muonium could be 
traveling at relativistic speeds, and it may be necessary 
to consider the effect of this boost on the wave functions 
[^ . Once observed, experiments focusing on precisely 
measuring the Lamb shift and hyperfine splitting (HFS) 
would put strong constraints on possible BSM solutions 
to the muon nroblem pH - l^ . 

Predictions for BSM models solving the muon prob¬ 
lem generically lead to corrections to the spectrum of 
true muonium as large as 0(100 MHz) (e.g, HI). This 
size corresponds to 0{ma^) corrections to the true muo¬ 
nium spectrum; therefore, we must first have a standard 
model prediction of this level. Currently, the theoreti¬ 
cal predictions for the HFS in true muonium are known 
fully to O(ma^) Additionally, all corrections of 

order 0{ma^) (see and references within for a his¬ 
torical review) and partial results for order 0{ma^) are 
known for positronium [stI - I^ . With an exchange of the 
electron mass for the muon mass, these results can be 
included in the true muonium prediction, to yield the 
current theoretical value of 

AUjil = 42330685(800) MHz, (1) 

where only the uncertainty arising from mo del-dependent 
hadronic effects [s^ are included. 

What remains of QED to be computed for true muo¬ 
nium at 0{ma^) are corrections that do not occur 
in positronium involving virtual electrons and hadrons. 
Without calculating all of these corrections, we would 
like some sense of the uncertainty, in the the¬ 

oretical value. To do so we make the following, rather 
gross, approximation. It is known that electron vacuum 
polarization corrections are the largest contribution from 
the work of unique to true muonium at 0{ma^). 

We will use this fact to estimate the unknown 0{ma^) 
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corrections unique to true muonium to be the complete 
0{ma^) diagrams multiplied by the photon polarization 
function, n(( 7 ^), that arises from the electron vacuum po¬ 
larization at momentum While not perfect, 

the logic behind this estimate is that insertion of electron 
loops into photon propagators of 0{ma^) diagrams con¬ 
stitute a large portion of the necessary corrections. This 
effect can estimated by the polarization function evalu¬ 
ated at the average expected momenta in the photon. For 
single annihilation photons, this means ~ pho¬ 
tons, and q^ ~ a^rrij^ for exchange photons. We will take 
the larger of these two. A final point in our estimate is 
that, since the majority of 0(ma^) contain two photon 
propagators, we should double the energy shift. 
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where Aif^g was obtained in and consists of the sum 
of the and terms in Table im Regardless of the 

accuracy of this estimate, these terms must be computed 
to obtain O(lOOMHz) precision. 

In this paper, we present a calculation of a novel con¬ 
tribution to the true muonium HFS that arises from the 
leading-order weak interactions, which scale as 0{ma^), 
and a more general expression for the leading pseudovec¬ 
tor contribution from a Z' particle, which have previously 
been considered as a BSM source in other atomic systems 

[Hil. 

Previous work on weak interactions in electronic sys¬ 
tems (ifil - lbll l have consistently found the effects to be 
heavily suppressed. This can be understood through the 
simple scaling associated with the weak interaction. For 
electroweak couplings in atomic systems, the dimension¬ 
less coupling constant is given by where ^ is the 

reduced mass of the system. In bound states, the 
term can be understood as arising partly from the wave 
function at the origin |'!/;(0)p = /irn^ where n is the 
energy level, and the remaining power of p, is needed for 
correct energy dimension. Taking the known value oi Gp, 
we see that the dimensionless constant can be expressed 
as Gp^i^ = 0.097fj?/M'^. From this result, an effective 
order of a can be computed for each atomic system via 

Gfm" = 0.097^ = a". (3) 


The leading-order electroweak correction to electronic 
systems should then be proportional to Gf|V'(0)P — 
Gpra^a^ ~ The theoretical uncertainties for 

the positronium, muonium, and hydrogen states are all 
far larger than this. 

This situation should be compared with that in true 
muonium and muonic hydrogen, where the coupling in¬ 
stead is proportional to Gpmf^ ~ leading to a naive 


prediction of its effect on the HFS of In muonic 

hydrogen, this effect will be difficult to discern because 
of the discrepancy in the proton radius and the unknown 
nuclear effects that plague hydrogenic atoms. In contrast, 
the almost purely QED atom of true muonium should 
offer a cleaner opportunity to measure these electroweak 
effects. As stated above, to see these effects, a number of 
QED contributions of 0{ma^) and 0{ma^) need to be 
computed. Putting these contributions together, a pre¬ 
diction for the HFS of true muonium to 1 MHz precision 
should be possible. 

In this work, we extend the previous work in [50l| by 
computing the leading-order energy shifts to HFS from 
the exchange and annihilation of a pseudovector in QED 
particle-antiparticle bound states, with special emphasis 
on true muonium, the most viable candidate in which 
to measure these effects in the near future. Numerical 
values using the Z boson are reported. In Sec. [TTl we re¬ 
view the previous calculations for the exchange process, 
and in Sec. uni we compute the results for the annihila¬ 
tion channel. We devote Sec. nYi to discussing the total 
contribution, and conclude in Sec. 0 


II. EXCHANGE OF Z BOSON 


Calculations have existed for many years for the ab¬ 
solute value of the contribution from the f-channel ex¬ 
change of a single Z boson to atomic systems [H, 113, S , 
but the correct sign for this effect was unclear. In 
the author resolved the issue, with the conclusion that 
the effect in hydrogen is positive, while in muonium it 
is negative due to the antiparticle nature of /r+. In [5l| 
the HES for other light muonic atoms were calculated, 
but the case of true muonium was conspicuously left un¬ 
touched. Eollowing [b^l, we write the leading-order La- 
grangian for the neutral weak interaction as 

Lint = gL^fiTzLZ^, (4) 


where g = gj cos0vy, and the weak left-handed spinor 
isodoublet L is defined using the convention: 

L = (5) 

While the Z boson has a vector component, this is fur¬ 
ther suppresed by a factor 1 — 4sin^ Ow Therefore we 
will neglect it’s contribution in this work. Generalizing 
the exchange correction result of for any particle- 
antiparticle bound state to the HES, we obtain 




2Gp {Zaf 3 

^/2 7rn® ^ 


( 6 ) 


where Z is the charge of the particle. Replacing the re¬ 
duced mass by mi/2 and taking Z = 1, we have that the 
correction to lepton-antilepton bound states is 


AA;“hfs(^) 


Gp mfa^ 

4-\/2 Trn® 


(7) 
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In a similar way, we may write the energy shift for a new 
Z' boson in the exchange channel as 




(n) 


mfa^ 
32M^, nn^ ’ 


( 8 ) 


where g is the new axial vector coupling between Z' and 
the leptons. 


III. ANNIHILATION VIA VIRTUAL Z BOSON 


For the true muonium and positronium systems, there 
is an additional diagram coming from the Z boson annihi¬ 
lation channel. Following , we find that the scattering 
amplitude for for this process is given by 

g^r3(2)r3(l) u(2)7^7^m(1)m(1)7^7^z;(2) 

4 (2m)2 - M2 ’ ^ 

where T^li) is the weak isospin of particle or antiparticle 
i. An additional minus sign occurs due to the relative sign 
in the weak neutral current of the antiparticle compared 
to the particle. This expression can be Fierz reordered 
to yield 

g2T3(2)T3(l) 

16[(2m)2 — M2] 

X [4z;(2)w(2)'u(1)m(1) -|- 2v{2)^^v(2)u{1)j^u{1) 

—2v(2)j^"f^v{2)u{l)jfj,j^u{l) — 411(2)7®t!(2)u(l)7®M(l)]. 

( 10 ) 


To simplify this expression, we take the non-relativistic 
limit (pi —^ 0), divide by the relativistic normalization 
'/2m for each spinor, and integrate over an s-state wave 
function to get the wave function squared at origin, mul¬ 
tiplied by constants. It is important to note that for 
pseudovectors with Mz> < mi , there are corrections that 
arise from the Yukawa-like interaction that depend on 
Mz'T where r is the radius of the atom. We will take 
the assignments of weak isospin for the antiparticle to be 
T 3 ( 2 ) = —1/2 and for the particle to be 13 ( 1 ) = 1/2. 
From this, we have the expression; 


( 11 ) 

Noting that (S”^) = we can rewrite the expec¬ 

tation value of the spins. Using the non-relativistic value 
for |'(/)(0)|2 given in Sec. U and evaluating this Hamil¬ 
tonian for the difference between the singlet and triplet 
state, (2 — S^) = 2, we have 


AKf. = 






8[(2to)2 — M|] TTO® 


( 12 ) 


For the case of a general Z' boson, this is the final ex¬ 
pression to consider after appropriate variable exchanges. 


It is interesting to note that for Mz' ~ {2m/, this term 
will dominate over the exchange term. This may poten¬ 
tially lead to stronger constraints on BSM pseudovectors 
for this mass range than those found for positronium 
in [ 2 ^. In contrast, for the real electroweak Z boson, 
we can safely approximate {2m/ — M^ k, —M^. Us¬ 
ing the definition that g from Sec. [TTl and the relations 
g2 = 47ra / sin^ , we have 


AKf, = - 


1 47ra 

8 sin^ 0w cos2 dwMj. 7rn® 


(13) 


Using the relation between the massive boson masses 
Mw = Mz cos{9w) gives 


1 47ra 
8 sin^ 9wM/ ttu® 


(14) 


Finally, using the relation Gf/v^ = 47ra/8M^ sin^(0u,), 
and that the lepton mass is related to the reduced mass 
by p = mi/2 we reach 


AKfe = - 


Gp mfa^ 
8 V 2 Trn® 


(15) 


Comparing this to Eq. 0, we see that the annihilation 
channel increases the electroweak correction by 50%. 


IV. TOTAL CONTRIBUTION 


Summing Eq. ([8|) and Eq. m, we find that the total 
pseudovector contribution to the hyperfine splitting of 
true muonium or positronium is 




+ 


4M|, 8[(2m)2 — M|,]J Sttu® ' 


mfa^ 


(16) 


Taking instead Eq. 0 and Eq. (US for the specific case 
of the electroweak Z boson: 


AEfhfs = - 


3Gf mfa^ 

8/2 TTn^ 


(17) 


Using our approximate value of Gpm/ ~ we find 
that as predicted, this correction does naively scale as 
0{m^/'). Closer inspection of the actual numerical 
coefficient reveals that this contribution is comparable 
in size to the 0{m^a^). For these particle-antiparticle 
bound states, the leading-order HFS is given by A.E^l = 
j^/mi. As a fraction of the leading order, we have cor¬ 
rections to true muonium of the size 

AU^,Ms = -2.58 X 10-®U® (18) 

and that the the total shift is —109 MHz. This correction 
has been included in Table m which gives the current 
state of the theoretical predictions. 
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TABLE II. Predictions for the Is HFS in true muonium. The designation of indicates corrections for true muonium missing 
from positronium, which depend upon m^/rrie. The error in the table consists of two estimated parts, the former from the 
model dependence of the hadronic loops, and the latter from electronic loop effects at 0{ma^), as described in text. * indicate 
terms that are only partially computed. 


V. CONCLUSION 

In conclusion, upcoming experiments present for the 
first time the strong possibility to detect and measure the 
properties of true muonium. These results will allow for 
strong discrimination of BSM explanations for the muon 
problem. To theoretically differentiate these models, one 
will need a precision of 0(100 MHz) for the hyperfine 
splitting of true muonium. 

In this paper, we have presented a subleading, but per¬ 
haps surprisingly large, correction from the electroweak 
interaction that needs to be accounted for in order to 
achieve the required precision. We have additionally 


computed the leading-order correction due to a general 
Z' boson. Further work remains to compute the more 
standard QED effects to higher order and to dramat¬ 
ically reduce the model-dependent uncertainties in the 
hadronic loops. 
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